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We propose a novel qnasiparticle interpretation of the eqnation of state of deconfined 
QGD at finite temperatnre. Using appropriate thermal masses, we introdnce a phe¬ 
nomenological parametrisation of the onset of conhnement in the vicinity of the phase 
transition. Lattice resnlts of bnlk thermodynamic qnantities are well reprodnced, the 
extension to small qnark chemical potential is also snccessfnh We then apply the model 
to dilepton prodnction and charm snppression in nltrarelativistic heavy-ion collisions. 

1. Introduction 

QGD is expected to nndergo a transition from a conhned hadronic phase to a deconhned 
partonic phase, the qnark-glnon plasma (QGP), at a temperatnre of Tc ~ 170 MeV [[^l- A 
central qnantity of matter in thermal eqnilibrinm is the Helmholtz free energy from which 
the pressnre p, energy density e and entropy density s - which are important ingredients 
for the description of nltrarelativistic heavy-ion collisions (URHIG) - can be derived. A 
first-principles nnderstanding of the eqnation of state (EOS) of hot QGD is therefore of 
great interest, also in order to reliably identify and calcnlate experimental signatnres of 
that elusive state. Perturbative results on the EOS are available up to order 0{g^) [|^, but 
show bad convergence for all temperatures of interest. Non-perturbative methods such as 
lattice QGD calculations hence become mandatory. From these numerical simulations the 
EOS of a pure gluon plasma is known to high accuracy [ and there are hrst estimates 
for the continuum EOS of systems including quarks [ii- 

Various interpretations of the lattice data have been attempted, most prominently as 
the EOS of a gas of quark and gluon quasiparticles. In a phenomenological framework, 
quarks and gluons are simply treated as non-interacting, massive quasiparticles [ ^]. More 
recently, a qnasiparticle description of QGD thermodynamics has been derived in a more 
rigorous treatment using resummed hard thermal loop (HTL) perturbation theory [ |^. 
Employing the full HTL spectral representions, the resulting EOS can be matched to 
lattice data down to temperatures T ~ 3 T^; below that temperature, non-perturbative 
physics not amenable in an expansion in Qs becomes important. Unfortunately, as evident 
from figure (left panel), current experiments only probe that very temperature regime 
where the underlying physics, the conhnement and chiral symmetry breaking mecha- 
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nism, is not sufficiently well understood. Phenomenological models incorporating as much 
physics as is known are therefore necessary. Here, we propose a new quasiparticle model 
of the QGP that incorporates a parametrisation of confinement close to T^, supplemented 
by thermal masses compatible with lattice results. Details can be found in [|^. 

2. Quasiparticles and confinement 

Consider a SU(3) gluon plasma. From asymptotic freedom, we expect that at very high 
temperatures the plasma consists of quasifree gluons. As long as the spectral function 
of the thermal excitations at lower temperatures resembles qualitatively this asymptotic 
form, a gluonic quasiparticle description is expected to be applicable. The dispersion 
equation of the gluonic quasiparticles reads 

ulc^e + mlin ( 1 ) 

where m^{T) acts as an effective mass generated dynamically by the interaction of the 
gluons with the heat bath background. 




Figure 1. Left panel; Sketch of theory vs. lattice and the temperatures probed by current 
experiments. Right panel: Normalised e, s and p (solid lines) in a quasiparticle approach 
[|^, compared to continuum extrapolated SU(3) lattice data (symbols) [||. 


However, the picture of a simple massive gas is presumably not appropriate close to 
Tc because the driving force of the transition, the conhnement process, is not taken into 
account. This physics has to be incorporated phenomenologically. Below Tc, the relevant 
degrees of freedom in pure SU(3) gauge theory are heavy, colour singlet glueballs. Ap¬ 
proaching Tc, deconfinement sets in and the gluons are liberated, followed by a sudden 
increase in entropy and energy density. Conversely, when approaching the phase tran¬ 
sition from above, the number of thermally active degrees of freedom is reduced due to 
the onset of conhnement. As T comes closer to Tc, an increasing number of gluons gets 
trapped in glueballs which disappear fiom the thermal spectrum: since mcB ~ 1-5 CeV 
and Tc ~ 270 MeV, glueballs are simply too heavy to become thermally excited in the 
temperature range under consideration (up to 5 Tc). So, all conhnement does statistically 
on a large scale is to cut down the number of thermally active gluons as the temperature 
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is lowered. This effect can be included in the quasiparticle picture by modifying the dis¬ 
tribution function of the gluons by a temperature-dependent confinement factor C(T): 



Figure 2. Left panel: The function B{T). Symbols show [H, the dashed line displays 
B{T) in other quasiparticle models [ Right panel: Normalised e, s and p (solid lines) 
for Nf = 2 -|- 1. The data points are from the lattice simulation in [ |^. 

To become quantitative, we have to specify the thermal masses m*(T) entering i). 
Based on the observation that the Debye screening mass tud evaluated on the lattice 
shows approximate critical behaviour [ Q , in accordance with a weakly first order phase 
transition and in contrast to perturbative results, we parametrise 

m.(T) ~ G„T (l - , (2) 

Consider now the entropy of a gas of massive gluons with such a dropping m^,{T). The 
result for s(T) will clearly overshoot the lattice entropy because light masses near 
lead to an increase in s(T). However, since the entropy is a measure for the number 
of active degrees of freedom, the difference may be accounted for by the aforementioned 
confinement process as it develops when the temperature is lowered towards T^. The 
explicit temperature dependence of the confinement factor C(T) can be obtained simply 
as the ratio of the lattice entropy and the entropy calculated with a dropping input gluon 
mass, and C{T) again shows near-critical behaviour: C{T) ~ (1 — Tc/T)'^. 
Thermodynamical quantities like the energy density can now be calculated as 

<T) = 16/^ [C(T)/bK)] uj, + B(T). (3) 

The function B{T) is not an independent quantity, but uniquely determined by m^{T), 
C{T) and their T-derivatives. It is necessary to maintain thermodynamical self-consistency 
and can be interpreted as the thermal energy density of the vacuum. The explicit ex¬ 
pressions for B{T), the pressure and the entropy density can be found in [ H. In figure |l| 
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(right panel), we compare e, s and p with lattice data, as a fnnction of T. We achieve a 
good and economic parametrisation. Fignre (left panel) shows the fnnction B{T) and 
the spacelike plaqnette expectation valne A^., as measnred on the lattice, that can be 
related to the thermal chromomagnetic condensate We hnd the simple relation 

B(T) = \a,{T)T‘ = 

with the zero-temperatnre condensate (G^)r=o- This correlation between B{T) and (B^)^ 
may hint at a deeper connection between B{T) as a carrier of non-pertnrbative effects, 
and the magnetic condensate. After all, B{T) represents the thermal energy of the (non¬ 
trivial) Yang-Mills vacnnm. 


3. Dynamical quarks 


The extension of the mechanism presented so far to systems with dynamical qnarks is 
not straightforward. Simnlations of fermions on the lattice are still plagned by problems. 
However, when plotting the lattice pressnre, normalised to the ideal gas valne, for the pnre 
gange system and for systems with 2, 2-1-1 and 3 qnark flavonrs, it is fonnd that the QCD 
EOS shows a remarkable flavonr independence when plotted against T/T^. The flavonr 
dependence is then well approximated by a term reminiscent of an ideal gas p{T, Nf) ~ 
(16 -|- 10.5 Nf) p{T/Tc) with a nniversal fnnctionp(T/Tc) [|T3- This hints at a conhnement 
mechanism being only weakly flavonr-dependent, and hence we assnme that the fnnction 
C{T) acts in a nniversal way on qnarks and glnons. Fignre ^ (right panel) shows pressnre, 
energy and entropy density for two massless qnark flavours and a heavy strange quark, 
compared to lattice data from a simulation with slightly different masses. The agreement 
especially close to Tc is certainly encouraging. 

For small quark chemical potential p, it is reasonable to assume that C{T) is only 
weakly p-dependent. The net quark density then takes the form 


ni{T,p) = 3j - fo] with / 


exp([a;fc =F/ii]/T) -hi’ 


( 4 ) 


where i labels the flavour. Figure ^ (left panel) shows Uq/T^ as a function of T/Tc, 
compared to preliminary continuum-extrapolated lattice data [ The general shape is 
well reproduced, which supports the validity of the quasiparticle assumption even close 
to Tc. Another quantity that tests the p-dependence is the quark number susceptibility 
x{T,pi)ij = dui/Opj. In hgure (right panel), we plot y = xiT,p = 0)^, normalised to 
the ideal gas value yo = NfT"^, along with corresponding lattice data [ 0 - Again, the 
agreement is satisfactory. 


4. Dilepton rates 

We apply the model now to dilepton production in URHIC where it enters in two ways: 
hrst, the dilepton emissivity of a static hot source, 

dMee ^ pi/(a;, fc;T) 
d^xdujd?k exp(a;/T) —l’ 


( 5 ) 
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Figure 3. Left panel: Net quark number density ng/T^ at small /i for 2+1 flavours. Data 
points are continuum estimates of lattice simulations with rriu^ — 65 MeV and ~ 135 
MeV [ Right panel: Quark number susceptibility for 2+1 flavours, normalised to the 
ideal gas value. The data points are continuum estimates of the lattice simulation of [ 

inn. 


is proportional to the vector spectral function py that is calculated from the photon self 
energy. By construction, the quasiparticle gg-loop is the only contribution in the QGP 
phase. In the hadronic phase, we use vector meson dominance to couple the photon to 
the = 1~ mesons p, uj and 0, taking into account temperature and baryon density 
effects 


12 


To compare with experiment, the rate (§) has to be convoluted with a hreball expansion. 
Assuming thermalisation, we use a setup reminiscent of hydrodynamics, constraining the 
hnal state by hadronic measurements and the initial state by geometry. The evolution 
inbetween is taken to proceed isentropically and fixed by the EOS of the quasiparticle 
model. We find an initial temperature of 300 MeV and a life time of the QGP phase of 7 
fm/c for SPS 30% central Pb(158 AGeV)+Au collisions. In the left panel of figure I, we 
show the agreement of the final result with the experimental data from the GERES/NA45 
collaboration. Data at a lower beam energy of 40 AGeV and pr-separated data are also 
well described |T2f . 


5. Charm suppression 


Next, we analyse J/il) suppression within kinetic theory [ |^], assuming that J/'ips 
are formed in the initial state of the collision and subsequently propagate through the 
thermalised medium. Interactions with the medium constituents lead to the break-up of 
the bound cc state. Neglecting coalescence, the time evolution of the J/V' number 
at midrapidity is then described by 


— N^ 


Y. ((^D^^rel))(r) pn{,T)N]j^ 


( 6 ) 


where Pn('r) is the density of the surrounding particles, and Urei))('r) an averaged 
cross section. Since the gluon density of the quasiparticle model is an order of magnitude 
larger than any hadronic densities, the dominant dissociation process is J/xjj + g —*• c + c, 
with no need to take into account hadronic comovers. Folding now the time-dependent 
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Figure 4. Left panel: Dilepton rate as a function of invariant mass Meg. Right panel: J/V’ 
suppression, normalised to the Drell-Yan yield, as a function of transverse energy Et- 

rate with the SPS hreball evolution (that is taken over unchanged from the dilepton 
discussion), we arrive at the right panel of hgure ^ Here, our result is compared with 
data from the NASO experiment. Again, the agreement with data is good in the region 
where thermalisation constitutes a valid concept. 

6. Summary 

We have constructed a quasiparticle model of the QGP near Tg with a phenomenological 
inclusion of conhnement, which works well for both gluons and quarks. Lattice data at 
small chemical potential could be naturally interpreted by the quasiparticle structures of 
our model. We then used the resulting realistic equation of state to construct a hreball 
and described both dilepton radiation and charm suppression within a unihed framework 
- a hrst step towards consistent heavy-ion phenomenology. 
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